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Abstract. The maximum force that a crawling cell can exert on a substrate is a quantity 
of interest in cell biomechanics. One way of quantifying this force is to allow the cell to 
^^ crawl against a measurable and adjustable restraining force until the cell is no longer 

able to move in a direction opposite to the applied force. Fukui et al. [Ij reported on 
an experiment where amoeboid cells were imaged while they crawled against an artificial 
gravity field created by a centrifuge. An unexpected observation was that the net applied 
1-^ force on the amoeba did not seem to be the primary factor that limited its ability to climb. 

l^^ Instead, it appeared that the amoeba stalled when it was no longer able to support a 

^Ny pseudopodium against the applied gravity field. The high g-load bend the pseudopodium 

T-H thereby preventing its attachment to the target point directly ahead of the cell. In this 

_-^ paper we further refine this idea by identifying the bending of the pseudopodium with the 

^~> onset of elastic instability of a beam under its own weight. It is shown that the principal 

^sj features of the experiment may be understood through this model and an estimate for the 

y—^ limiting g-load in reasonable accord with the experimental measurements is recovered. 
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Motility is a fundamental trait that distinguishes living things. The macroscopic motion 
of plants and animals can ultimately be reduced to motion on the level of single cells |2] . 
The movement of single cells also play a pivotal role in phenomena not directly related to 
movement of the whole organism. Thus, the morphogenetic migration of cells is responsible 
for the appearance of form and structure in embryogenesis. The movement of cells is 
crucial in the process of wound healing and in the functioning of the immune system and 
unfortunately in metastasis, where cancer cells spread from the primary tumor to invade 
other organs of the body |i3j. When the environment of the cell is a fluid, the cell swims, but 
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on a substrate or in a fibrous network, cells crawl. The mechanics of crawling [H O [6l [71 18] 
is a fundamental problem in biomechanics that has not yet been fully understood. 

The measurement of the actual forces exerted by crawling cells on substrates is of obvious 
importance and various techniques have been employed in order to perform such mea- 
surements [9]. One method is to restrain the cell in some way, for example by holding 
a micro-needle in its path or by applying a suction force with a micro-pippette. An in- 
genious non-invasive technique involves imaging the wrinkles on an elastic substrate from 
which the applied force may be inferred [TU]. Restraining forces can be applied to crawling 
cells by attaching magnetic beads to them and pulling with magnetic fields. Gravity is an 
excellent candidate for an external force as it does not require anything to be attached to 
the cell and it is the least likely to affect the natural behavior of the cell. Fukui et al. [1] 
reported on an experiment in which amoebae of Dictyostelium discoideum were allowed to 
crawl against an artificial gravity field created by a centrifuge. They determined the max- 
imum g- force at which the amoeba "stalled" - that is, was unable to crawl in a direction 
opposite to the gravity field. However, a rather surprising finding was that the net force 
alone did not limit the crawling ability of the amoeba. Indeed, when the culture medium 
was replaced by one of a higher density, so that the amoeba was actually buoyant in it, 
it nevertheless stalled, even though, in this case, the gravity field was pushing the cell in 
the direction of motion. Fukui et al. observed that what appeared to limit the crawling 
ability of the amoeba was not the net force on the cell but rather the inability of the cell to 
extend a pseudopodium against the large gravitational field generated by the centrifuge; 
the pseudopodium was observed to bend and therefore not able to attach to a point on the 
substrate directly ahead of the amoeba (see Figure 1). In the present study, we show that 
this proposal of Fukui et al. that the buckling instability of the pseudopodium is the major 
factor limiting the ability of cells to crawl against high gravitational forces, is supported by 
a quantitative analysis based on mechanics. We further suggest that the ability to crawl 
against high gravitational forces is significantly improved in cells that are able to actively 
generate an internal turgor pressure in the pseudopodium by virtue of the contractile forces 
in the posterior part of the cell cortex [7J . 

If one assumes, based on these observations, that the structural stability of the actin 
rich pseudopodium against the gravity field limits the crawling ability of the cell, then 
it should be possible to use classical results on the stability of elastic structures under 
gravity to estimate the stalling acceleration {gc)- The relevant result is the classical Euler- 
Bernoulli theory which predicts that an elastic beam of uniform cross-section buckles under 
a compressive force [F) if this force exceeds a critical value given by 

(1) F=^ 

where E is the Young's modulus of the material, I is the area moment of inertia of the 
cross-section, L is its length and the value of the numerical constant a depends on the 
conditions at the ends of the beam; for a beam clamped at one end and free at the other, 
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Figure 1. The amoeba stalls (is unable to move "up" against the artificial 
gravity field g of the centrifuge) when the pseudopodium advancing in the 
direction of intended motion is bent by gravitational forces. Thus, the 
pseudopodium fails to attach at the target location A on the substrate 
directly above the cell body. Downward or lateral crawling of the cell is not 
prevented. 



a = 7r^/4. The Euler-Bernoulli theory proceeds from the assumptions of mechanical equi- 
librium and small deformations that result in a linear boundary value problem for the 
beam center line. The requirement that this equation admit nonzero solutions results in an 
eigenvalue problem, and Eq. [II corresponds to the lowest eigenvalue [llj. If the compres- 
sive force is the weight of the beam (W), a simple estimate for the maximum weight up 
to which the vertical configuration could be stable may be obtained by assuming all of the 
weight to be concentrated at the center of mass, so that 

(2) H' = ^. 

A more careful analysis where the weight is assumed to be uniformly distributed along the 
beam was provided by Greenhill [12]. Eq. [2] was found to hold except that a ^ 1.99. For 
tapered beams Eq. [2] may still be applied if I is regarded as the area moment of inertia of 
the base. Greenhill showed that for a right circular cone a ~ 2.54 and for a paraboloid of 
revolution a ~ 2.47. Keller and Niordson [13] calculated the greatest height L that a beam 
of fixed weight W can have if it is allowed to taper in an arbitrary way while preserving 
the shape of the cross-section. Once again, the result can be expressed in the form of Eq. [2] 
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with a ~ 8.23. Thus, Eq. [2j where a is a numerical constant that has a value roughly 
between one and ten, may be used to estimate the critical acceleration Qc of the centrifuge. 
For this purpose, one can rewrite Eq. [2] as 

4aEI 

using the apparent weight [16j (in place of W) of the pseudopodium which is the difference 
in density between the pseudopodium and the culture medium (Ap) multiplied by the 
volume (V) of the pseudopodiurrr) 

If the pseudopodium is regarded as a circular cylinder of diameter 2 fim and length 5 /im, 
then / ss 0.8 X 10~^^ m^ and 1/ ?a 1.6 x 10~^^ m^. Fukui et al [I] estimate that the density 
of the actin rich pseudopodium must be at least p = 1.124 gm/cm^. If we accept this value, 
then, since the medium density (at % Percoll) is po = 1.005 gm/cm^, Ap = p — po = 0.119 
gm/cm^. The greatest uncertainty arises in estimating the Young's modulus E. If the 
pseudopodium is presumed to be supported predominantly by the mechanical strength of 
the actin network one could use in vitro measurements of the strength of actin gels. Such 
measurements are usually expressed in terms of the shear modulus G which is related to 
E and the bulk modulus Ky as E = 9GKy /{2>Kv + G) ss 3G (since G/Ky is smaller 
than 10"'''). In vitro measurements show that G for actin networks vary from about 300 
Pa in the "gel" state to a value three orders of magnitude lower |14| [T5] in the fluid state. 
The measured value depends primarily on the density and length of actin filaments and 
the density of cross- links created by various actin binding proteins (ABPs). If we take 
G w 300 Pa and a w 2.5 (corresponding to a structure of paraboloid shape) Eq.[3] yields a 
numerical estimate (7c ~ 5 x 10^ m/s^ . 

In the experiment pQ Fukui et al. found that the myosin II knockout mutant (HSl) of the 
Dictyostelium amoeba stalled at Qc ~ 3.9 x 10^ m/s^ (in the buffer without Percoll) which 
is in reasonable accord with the above estimate. The wild type strain (NC4) containing 
myosin II did not stall even at the highest accelerations tested (about 11.2 x 10^ m/s^). 
This is probably because the wild type cells are able to create considerable turgor pressure 
due to the myosin Il-dependent contractile forces in the actin cortex thereby stiffening the 
pseudopodium. This case is discussed next. 

In order to understand the effect of turgor pressure one needs to recognize that the pseu- 
dopodium is a poroelastic solid. When a beam made of such a material bends, the compres- 
sion of the pores on one side of the neutral line results in a pressure rise in the interstitial 
fluid which then drives a flow across the beam. On account of viscous resistance, the fluid 



this correction in the formula for the stability threshold is necessary because the equation for the defor- 
mation y{x) of the beam is obtained by minimizing the energy functional L dx [EI{y")^ /2 + (p — po)gSX^ 
where 5X is the vertical displacement of a point on the beam at location 'x' due to beam curvature and 
Po is the density of the external medium. This differs from the corresponding expression in the absence of 
the external medium in that the density difference Ap — p — po replaces the density of the material of the 
beam, p. 
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Figure 2. Diagram illustrating the radial expansion of the pseudopodium 
as a consequence of the rise in internal turgor pressure, resulting in a "stiff- 
ening" of the structure against elastic buckling. 



responds to the bend with a time lag so that the system behaves much like a coupled mass- 
spring-damper system [T7j. However, this mechanism does not alter the stability limit since 
the bifurcation at the onset of elastic instability takes place at zero frequency. Internal 
hydrostatic pressure stiffens the structure by one of two mechanisms (a) the swelling and 
consequent stretching of elastic elements may put it in a regime where the stress strain 
relation is no longer linear (b) the swelling may alter geometrical parameters (specifically, 
the parameter / in Eq. ^. The first of these effects is less likely, though a careful esti- 
mate is difficult as not much information is available on the nonlinear elasticity of actin 
gels. However, the second of these effects is readily estimated. The situation is depicted 
in Figure 2. Since an element of the cell membrane is in equilibrium due to the balance of 
an outward pressure (Ap) and an inward elastic stress E{6R/R) where 5R is the increase 
in the radius R, we have 5R/R = {Ap)/E. Thus, 



(4) 



1 + 



9c 



1 + 



5R 
~R 
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since the area moment of inertia of a cylinder is / = irR'^/A. In the experiment [T] the 
wild type strain (NC4) containing myosin II did not stall at the maximum acceleration of 
gc = 11.2 X 10"^ m/s^, suggesting that 6gc/gc > 1-9. Thus, Ap > 0.305 x E = 275 Pa, 
using the value E = 3G = 900 Pa cited earlier. Pasternak et al. [18j report a difference 
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AT w 0.13 mdyn per micron in the cortical tension between the strains AX4 and the 
myosin lacking mhcA- strain of Dictyostelium both in the resting phase. If one converts 
this number to a pressure using the Laplace formula for surface tension, one obtains the 
estimate for the myosin generated pressure: Ap = 2AT/R ~ 520 Pa where i? ~ 5 /xm is 
taken as a characteristic radius of the cell. Thus, our rough estimate Ap > 275 Pa is not 
inconsistent with reported values for myosin Il-dependent pressures in the cell cortex that 
may be inferred on the basis of existing experimental data. 

In conclusion, the hypothesis advanced by Fukui et al. that the ability of the Dictyostelium 
amoeba to support a pseudopodium against a strong gravity field limits its ability to crawl 
against such a field appear to be supported by the present analysis based on the mechanics 
of the buckling of elastic structures. Similar ideas have long been used in the field of plant 
biomechanics [19j, but the centrifuge experiments of Fukui et al. present an opportunity 
for an application of these concepts to the mechanics of cell crawling. In a broader context, 
we illustrate that the interpretation of experiments in which one attempts to quantify the 
force applied by a crawling cell on a substrate by measuring an applied restraining force |9] 
may be subtle, because, the limiting factor may not be the ability of the cell to pull against 
the applied force but rather a failure in some other aspect of the motility mechanism of 
the cell (e.g. contact inhibition |20] ) . 
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